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Abstract 

This paper studies the Vlasov-Monge- Ampere system (VMA), a 
fully non-linear version of the Vlasov-Poisson system {VP) where the 
(real) Monge- Ampere equation det gf-^r ■ = p substitutes for the usual 
Poisson equation. This system can be derived as a geometric approx- 
imation of the Euler equations of incompressible fluid mechanics in 
the spirit of Arnold and Ebin. Global existence of weak solutions and 
local existence of smooth solutions are obtained. Links between the 
VMA system, the VP system and the Euler equations are established 
through rigorous asymptotic analysis. 

1 Introduction 

The classical Vlasov-Poisson [VP) system describes the evolution of an elec- 
tronic cloud in a neutralizing uniform background through the following equa- 
tions 

(1) l^ + e- v,./ + v,<^- V5/ = 

(2) e^A^ = p - 1, 

where f(t, x,^) > denotes the electronic density at time t > 0, point x G M*^, 
velocity ^ G M'^ (usually d = 3), p(t,x) > denotes the 'macroscopic' density 

(3) p(t,x)= / /(t,x,Oc?e, 

and (p{t,x) denotes the electric potential at time t and point x generated, 
through the Poisson equation (j2|), where e is a coupling constant, by the 
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difference between the electronic density p{t, x) and the neutrahzing back- 
ground density, which is supposed to be uniform and normahzed to unity. 
Standard notations V = ■■■,dd) and A = df + ... + d'^ have been used 
and ■ stands for the inner product in R''. The mathematical theory of the 
VP system is now well understood. In particular, existence of global smooth 
solutions in three space dimensions has been proved in (2^ (see also jTH]; 
|26j). In the present paper, a fully nonlinear version of the VP system is 
addressed : 

(4) ^ + e-V./ + V.</^- V€/ = 

(5) det{I + e^D^ip) = p, 

where the (real) Monge- Ampere equation substitutes for the Poisson equa- 
tion (j2)). Here, D'^(p{t,x) stands for the d x d symmetric matrix made of all 
second order x— partial derivatives of ip, I stands for the d x d identity ma- 
trix and det for the determinant of a square matrix. The occurrence of the 
Monge- Ampere equation in mathematical modeling is not very common. No- 
tice, however, that a very similar system can be found in meteorology with 
Hoskins' semi-geostrophic equations (cf. [HI, [IB] and the included references). 
In a simplified two dimensional setting, the semi-geostrophic equations read 

(6) |^ + Mp} = 

(7) det{l + e^D^ip) = p, 
where {■, ■} denotes the usual Poisson bracket. 

Formally, as the coupling constant e is small, the VP and VMA equations 
asymptotically approach each other up to order 0{e^). Indeed, linearizing 
the determinant about the identity matrix leads to 

(8) det{I + e^D^^) = l + e'^Aip + 0{e*). 
The formal limit, as e = 0, reads 

(9) ^ + V./ + V.</^- V5/ = 

(10) p=l, 

where constraint ()1U|) substitutes for both the Poisson and the Monge- Ampere 
equations. The limit system ()9|10j) . that we call constrained Vlasov system, 
can be seen as a 'kinetic' extension of the Euler equations of classical incom- 
pressible fluid mechanics, 

(11) dtv + {v ■V)v = -Vp 

(12) V-v = 0, 
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where v{t,x) G M'^ and p{t,x) G M respectively are the velocity and the 
pressure of the fluid at time t and position x. Indeed, any smooth solution 
{v,p) provides a 'monokinetic' solution to the constrained Vlasov system 
dnCni), defined by 

f{t,x,^) = 5{^-v{t,x)), (p = -p. 

Here a monokinetic solution means a delta-valued solution in the ^ variable. 
In addition, the constrained Vlasov system ()9I10|) turns out to be a natural 
extension (or F limit) of the Euler equations from both geometrical and 
variational reasons, as explained in section |2l 

In a similar way, there is a monokinetic version of the VP system, the so- 
called (pressureless) Euler-Poisson (EP) system, which reads 

(13) dtv + {v ■ V)v = Vif 

(14) dtp + V- (pv) = 

(15) e'^Aip = p-l. 

A rigorous asymptotic analysis of the VMA system as e — > will be provided 
f sections 15.11 and 15 .2^ . in the case when the initial electronic density 

(16) /(t = o,x,o = r(x,o 

is asymptotically monokinetic, namely approaching — vq{x)), for some 
smooth divergence free velocity field f o, as e tends to zero. Before this asymp- 
totic analysis, we want to explain the geometric origin of the VMA system. 
It has been known, since Arnold's celebrated work (cf. [2 ), that the Eu- 
ler equations (formally) describe geodesies curves along a suitable group of 
volume preserving maps, lengths being measured in the sense. We will 
show (section 121) that the VMA system just describes approximate geodesies 
obtained through a very natural penalty method, where e stands for the 
penalty parameter. For this geometric interpretation to be valid, the Monge- 
Ampere equation (jSJ must be understood in the following weak sense: for 
each fixed t, (f{t, ■) is the unique (up to an additive constant) function such 
that \E'(a;) = x^/2 + e^v9(t, x) is convex in x and 

(17) \/geC%R^),[ g{V^{x))p{t,x)dx = [ g{y)dy, 

where f2 is a fixed bounded open convex set where the neutralizing back- 
ground of the electrons is assumed to be located. (This definition is made 
precise in section EISl) Notice that, by construction, V\E' must be valued in 
the closure of fl and, therefore, the potential (p enjoys the following property 

\x + e^Vxfit, x)\ < sup \y\ < +oo. 

yen 
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There is no similar bound for the electrostatic potential of the classical VP 
system. Thus, in some sense, the VMA system can be seen as a nonlinearly 
saturated version of the VP system. 

Beyond the geometric derivation of the VMA system, our main analytic re- 
sults follows: 

• The VMA system admits global energy preserving weak solutions. 

• The VMA system admits local strong solutions in periodic domains. 

• For well prepared, nearly monokinetic initial data, the solutions of the 
VMA system converge when e goes to to those of the Euler equations. 

• In this asymptotic, the EP system is a higher order approximation of 
the VMA system. 

The paper is organized as follows: in section 2, we first recall the geometric 
nature of the Euler equations, then we explain why the constrained Vlasov 
system fl9lTn|l is a natural extension of the Euler equations from a variational 
point of view, finally we introduce the concept of approximate geodesies for 
volume preserving maps, and derive the VMA system. Section 3 is devoted to 
the proof of existence of global energy preserving weak solutions. In section 4, 
we prove existence of local strong solutions, in the case of a periodic domain. 
Finally, in section 5, we study the asymptotic behavior of the VMA system 
as e goes to 0. 

2 The geometric origin of the Vlasov- Monge- 
Ampere system 

2.1 The Euler equations 

The motion of an incompressible fluid in a domain C M'^ is classically 
described by the Euler equations (E): 

(18) dtv + {v ■ V)v = -Vp 

(19) V-w = 0, 

with t & W, X & Q, where v = v{t,x) stands for the velocity field and p = 
p{t, x) for the scalar pressure field. These equations have a nice geometrical 
interpretation going back to Arnold (see p|). Introducing G{Q) the group 
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of all volume preserving diffeomorphisms of with jacobian determinant 
equal to 1, and measuring lengths in the sense, we may define (at least 
formally) geodesic curves along G{Q). It turns out that the Euler equations 
just describe these curves. For the same reasons, the Euler equations can be 
seen as the optimality equations for the corresponding minimization problem: 
given two maps chosen in G{Q), find an L^— shortest path between them 
along G{fl). It was shown by Shnirelman (see also |2] and ^8]) that, 
in the case when Q is the unit cube in M^, there are many maps for which 
there are no such shortest paths. Beyond this negative result, [Hj established 
that minimizing paths are more appropriately described by doubly stochastic 
measures. These measures (also called polymorphisms) generalize volume 
preserving maps in the following way: a doubly stochastic measure fi{dx, dy) 
is a (Borel) probability measure on 1] x with two projections on each copy 
of Vt both equal to the (normalized) Lebesgue measure. It is known -see |22] . 
for instance- that any such /i can be weakly approximated by a sequence 
fj,n{dx, dy) = 6{x — gn{y))dy where each is a volume preserving map of fl. 
In inj it was shown that, in the case considered by Shnirelman for which there 
is no classical shortest path, minimizing paths along G{Q) converge to paths 
of doubly stochastic measures t — > fi{t; dx, dy) governed by the following 
extension of the Euler equations 



where v = v{t; x,y) and p = p{t,x) can be respectively seen as the velocity 
field and the pressure field attached to /i. (Notice that the velocity field v 
generally depends on the extra variable y and is not a classical but rather a 
multivalued velocity field.) These equations are just a reformulation of the 
constrained Vlasov system ()9ll0p . Indeed, it can be checked, under appro- 
priate regularity assumptions, that the kinetic measure / defined by 



solves ()9I10|) when {fi,v,p) solves ()20I21|) . Thus we conclude that the con- 
strained Vlasov system ()9|10|1 is a natural variational extension of the Euler 
equations. 

2.2 Approximate geodesies 

A general strategy to define approximate geodesies along a manifold M (in 
our case M = G{Q)) embedded in a Hilbert space H (here H = L^(f2,M°')) 



(20) 
(21) 



dtfi + ■ (yUf ) = 0, 




(22) 
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is to introduce a penalty parameter e > and the following unconstrained 
dynamical system in H 



(23) duX + ^Vx {d\X, M))) = 0. 

In this equation, the unknown t — > X{t) is a curve in H, d{X,M) is the 
distance (in H) of X to the manifold M, i.e. in our case as M = G{Q), 

(24) d{X,G{Q))= inf \\X - g\\H, 

and, finally, Vx denotes the gradient operator in H. This penalty approach 
has been used for the Euler equations by the first author in [Zj. It is similar- 
but not identical- to Ebin's slightly compressible flow theory ^Hj, and is a 
natural extension of the theory of constrained finite dimensional mechani- 
cal systems The penalized system is formally hamiltonian in variables 
{X, dtX) with Hamiltonian (or energy) given by: 

E = \\\dtx\\i + ^^d\xMm 

(Multiplying equation (PHj) by dtX, we formally get that the energy is con- 
served.) Therefore it is plausible that the map X{t) will remain close to G{VL) 
if properly initialized at t = 0. A formal computation shows that, given a 
point X for which there is a unique closest point Tlx to X in the H closure 
of G{yt), we have: 

(25) Vx (diX, G)) = ^^^(^ - vrx). 
Thus the equation (j^Hj) formally becomes: 

(26) dttX + \{X - nx) = 0. 

To understand why solutions to such a system may approach geodesies along 
G{Q) as e goes to 0, just recall that, in the simple framework of a surface 5* 
embedded in the 3 dimensional Euclidean space, a geodesic t s(t) along S 
is characterized by the fact that for every t, the plane defined by {s(t), s{t)} 
is orthogonal to 5". In our case, dttX(t) is nearly orthogonal to G{fl) thanks 
to ()26|). meanwhile X{t) remains close to G{fl). 

The approximate geodesic equation was introduced in in order to allow a 
spatial approximation of G{fl) by the group of permutations of N points Aj 
chosen to form a discrete grid on Q. On such a discrete group, the concept of 
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geodesies becomes unclear meanwhile approximate geodesies still make sense. 
They can be interpreted as trajectories of a cloud of N particles Xj moving 
in the Euclidean space M'^^, which substitutes for H. These particles solve 
the following coupled system of harmonic oscillators 

where a is a time dependent permutation minimizing, at each fixed time t, 
E — ^cr(i) among all other permutations of the first N integers. The 
convergence of this discrete model to the incompressible Euler equations for 
well prepared initial data was proved in [Zj . In order to study the continuous 
version (|26|l . a specific study of the projection problem (|2l| is needed. 



2.3 The polar decomposition Theorem 

Let us first recall a general measure theoretic definition: 

Definition 2.1 Let A and B he two topological spaces, let p he a Borel finite 
measure of A and X a Borel map A ^ B , we call the push-forward of p hy 
X and note X^dp the Borel measure rj on B defined hy 

V/ G C\B), [ f{y)d^{y) = [ f{X{x))dp{x). 
Jb J a 

Let us now consider the case of a bounded open subset fl of the Euclidean 
space MJ^ equipped with the Lebesgue measure that we denote dx. We say 
that a Borel map s : — > f2 is volume (or Lebesgue measure) preserving if 
s^dx = dx, i.e. if for all g G C°(fi) one has J^g{x)dx = J^g{s{x))dx, or 
equivalently, for any Borel subset B of Q one has |s~^(S)| = \B\. The set 
of all measure preserving maps of is a closed subset of the Hilbert space 
H = L'^{Q, W^) and will be denoted by S{fl). Notice that S{Q) is only a semi- 
group for the composition rule and contains the group of volume preserving 
diffeomorphisms G{Q). It is known [23j that, at least in the case when Q is 
convex and d>2, S{n) is exactly the closure of G{n) in L\n,M.'^), which 
implies d{.,G{n)) =d{.,S{n)). 

The polar decomposition Theorem for maps jH] (extended to Riemannian 
manifolds in [21]) will be crucial for our analysis of the VMA system: 

Theorem 2.2 Let^l he a hounded convex open suhset ofMJ^, letX G L^(fi; W^) 
and px = X^dx, where dx is the Lehesgue measure on fl. Assume px to he a 
Lehesgue integrahle function, or, equivalently, X to satisfy the non- degeneracy 
condition: 

(27) ME C Borel , |^| = ^ \X''^{E)\ = 0. 
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Then there exists a unique pair {W^XiT^x) where $x o convex function 
and nx € S{fl), such that 

(28) X = V$x o vTx. 

/n this 'polar decomposition', Tlx is also characterized as the unique closest 
point to X on S{Q) in the sense and $x is characterized to be (up to an 
additive constant) the unique convex function on Q satisfying 

(29) / 9{x)dpx= I g{X{y))dy= [ g{V<l>x{y))dy, 

for any g G CO(M^) such that \g{x)\ < C{1 + \x\'^). 

In addition, the Legendre-Fenchel transform of $x defined by 

(30) ^x (x) = snp{x-y- ^x (y) } 

yen 

is Lipschitz continuous on W^, with Lipschitz constant bounded by sup^gj^ |x| 
and has the following properties : 
V'^xix) G fl holds true for px a.e. x, 

(31) / giV^x)px{x)dx = [ giV^x{X{x)))dx = I g{x)dx 



JK.<i Jn Ju 

for any g G C^{Q), and 

(32) V^xC^'^xix)) = X px{x)dx a.e, 

(33) V^x(V$x(l/)) =?/ dya.e, 

(34) 7rxiy)=V^xiX{y)) dya.e. 



We make here several remarks on Theorem 12.21 

Link with the Monge-Ampere equation We can interpret (j29|l as a weak 
version of the Monge-Ampere equation: 

px(V$) detD^<l> = 1 

and ()3ip can be seen as a weak version of another Monge-Ampere equation: 

det = Px 

V\l/ maps supp(px) in fi. 

The pair ($X) ^x) depends in fact only of Q and the measure px = Xjj^dx, 
and if condition fails, then existence and uniqueness of the projection 
Tlx may fail, but existence and uniqueness of V<I>x remain true. 

Theorem 12.21 and the subsequent remarks allow us to introduce the fol- 
lowing notation that will be used throughout the paper: 
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Definition 2.3 Let Q be a fixed bounded convex open set of M.'^, let p be 
a positive measure on of total mass absolutely continuous w.r.t the 
Lebesgue measure and such that /(I + \x\'^)dp{x) < +00. We call p], or, 
in short, $[p], the unique up to a constant convex function on Q satisfying 

(35) WgeC\m'')nL\dp), [ g{x)dp{x) = [ giV^Q, p]{y))dy. 

We call \l/[fi,p] its Legendre-Fenchel transform satisfying 

i36)\/geC\R')nL\n,dy)J giV^n, p]ix))dpix) = !g{y)dy. 

If no confusion is possible we may write $ (resp. instead of p] (resp. 

We will use some additional results from The first one establishes the 
continuity of the polar decomposition: 

Theorem 2.4 Let p be a Lebesgue integrable positive measure on M."^, with 
total mass fl, such that /(I + \x\'^)dp < +00. Let pn be a sequence of 
Lebesgue integrable positive measures on W'' , with total mass Vt, such that 
Vn, /(I + \x\'^)dpn < +00. Let = $[fi,p„] and \E'„ = ^[Vt, pn] be as in 
Definition\2~B If for any f G CO(R'^) such that \f{x)\ < C(l + Ixp), / / rfp„ 
converges to j f dp , then 

• $„, converges to p] uniformly on each compact set ofQ and strongly 

• \E'„ converges to \E'[fi, p] uniformly on each compact set o/M'^ and strongly 
in W^'^{K) for every K compact in W^. 

The second one provides a 'dual' definition of the distance between a map X 
and the semi-group SiVt): 

Theorem 2.5 Let X G R'^) and p = X^dx, where dx is the Lebesgue 

measure on Q. Assume p to be a Lebesgue integrable function. Then 

^d\X,Sm = J {\x\'/2~^n,p]{x))p{x)dx + Jj\y\y2-^n,p]{y))dy 

= sup / {\x\'^/2 ~ u{x)) p{x)dx + / {\y\'^/2-v{y))dy, 
u,v J Jq 

where the supremum if performed over all pairs {u, v) of continuous functions 
on W'- such that u{x) + v{y) > x ■ y pointwise. 
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2.4 The Vlasov-Monge- Ampere system 

Let us now derive the VMA system as the kinetic formulation of the approx- 
imate geodesic equation (j26|l . First, from the polar decomposition Theorem 
\2.2\ equation (j26|l reads 



(37) duX{t,x) = V^{t,X{t,x)), 

where 

vnn,p{t,-)]{x)-x 



(38) ^^{t,x) 



e2 



and \E'[i7,p] is as in Definition ()2.3|) . This means that satisfies (0) in a 
weak form with the additional condition that the range of x — *• x + e'^'Vip{t, x) 
is contained in Q. 

Next, let /° > be a given initial density function, that we assume to 
be in L°°{M.'^ x Mf^), compactly supported and satisfying the compatibility 
condition 



(39) J f{x,Odxd^=\n\. 

For each t > 0, let us define (x,^) f{t,x,^) to be /° pushed forward by 
the following ODE 

(40) dtX{t,x,0 = ^it,x,0 

(41) d,E{t,x,0 = (yy^)iX{t,x,0) 

(42) (X,S)(t = 0,x,O = (x,O- 

Then / satisfies the following kinetic (or Liouville) equation 

(43) % + + ^« ■ = ° 

(44) /(0,.,.) = r, 

which must be understood in the following weak sense 

I / ( + ^ V.(? + Vy^ ■ V^g] fdxd^ 

(45) =-[ fo{x,09{t = 0,x,0dxd^. 
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This linear Liouville equation is nonlinearly coupled to equation (jHHj) . where 
p is hnked to / by equation Q. Finally, we have defined, through ()38|43|44|) . 
the weak formulation of the VMA initial value problem. 
The energy of the system is defined by 



3 Existence of global renormalized weak so- 
lutions 

The main result of this section is as follows: 

Theorem 3.1 Let ^ P{x,i) > be m L°°(R'^ x R''), with compact 

support in both x and ^, satisfying condition K^) . 

Then the VMA system \38\4^44^ admits a global weak solution (/, p, \E'), 
with f G X X R'^) and (p, Vtp) G L°°([0, T] x R'^) for allT>0. In 

addition, each such weak solution enjoys the following properties: 

• f is a continuous function oft, valued in L^(R'^ x R'^), for every 1 < 
p < oo 

• the density p is a continuous function oft, valued in L^(]R'^), for every 
I < p < oo, 

• the support of f{t, -, ■) in (x, ^) is compact, with a diameter growing no 
more than linearly in t. 

• the total energy defined by is conserved, 

• the 'renormalization' property (in the sense of mV 



(46) 




dg{f) 



+ V.-(e^7(/)) + V5-(VvP^7(/)) = 



dt 



holds true for all g G C 



(M) 



the trajectories of \4^4^ ^'^c uniquely defined for almost every initial 
condition {x,C,), 



t — > /(t, ■, ■) is just pushed forward along the trajectories of (^TjjT^ 
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Proof of Theorem 13. It 

We build a sequence of approximate solutions {fh, ^h)h>o by time discretiza- 
tion and let the time step h go to zero. To handle the limiting process, the 
non-linear terms will be treated with the help of Theorem 12 .41 More precisely 
if one can extract a subsequence such that, for every t, fh(t, •, •) converges 
weakly, then we can deduce from Theorem 12.41 that the corresponding se- 
quence V\l//i(t, ■) will converge strongly, and this will allow us pass to the 
limit in the nonlinear term. 



3.1 Construction of a sequence of approximate solu- 
tions 

We consider rj G C^{W'-) such that rj > 0, J^^V = ^ Vh = -^vih)- We 
then seek approximate solutions as solutions of the approximate problem 

(47) f + e-v./. + ^M^.v,A = o 

(48) /,(0, X, = A°(x, = /o Vh ® Vh 

(49) "^hit) =Vh* ^[fi, pit = nh)] for t G [nh, {n + l)h[. 

V^/i being a smooth function of space this regularized equation admits a 
unique solution that one builds by the method of characteristics. Since the 
flow is divergence-free in the phase space, the solution satisfles 

(50) Vp G [l,+00], ||/h,(t)||LP(IRdxIRd) = ||A(0)||LP(RdxIR<')- 

By construction (through Theorem l2.2|l . V\E'/i is valued in the convex bounded 
set fl. Suppose that /°(x, ^) vanishes outside of the set {x^ + e^^^ < C^} for 
some constant C > flxed and denote R = sup^^j^ \y\. Then we have 

Lemma 3.2 Vt > 0, fhit, -, ■) zs supported zn {^/x^T^ <C + Rt/e}. 
Proof : We just write 

e^duX + X = V^j,{X) 

in complex notation —iedtZ + Z = F, where Z = X + iedtX and F = 
V^h{X), which is bounded by R. This leads to 

Z{t) = Z{0) exp{-it/e) + ie~^ [ exp(-i(t - s) /e)F{s)ds 

Jo 

ant the desired bound easily follows. Notice here a sharp contrast with the 
classical VP system, for which the ^—support of the solutions cannot be 
controlled so easily (except in the one dimensional case). □ 
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Convergence of the sequence of approximate solutions 
Using ()50|) and Lemma 13.21 there exists, for any 1 < p < cxd, up to the 
extraction of a subsequence, / G L^([0,T] x M'^ x R"') such that fh converges 
weakly to / as /i ^ 0. 

It remains to show that the product fh^^h converges to the good limit. 
For this we need strong convergence of V'^h- We already know that V\l//i G 
L°°([0,T] X R'^). We claim that for all t > 0, V^hit,-) converges strongly 
to V\E'(t, ■) in Ll^^iW^), Wq G [1, +oo[. Indeed, such a strong convergence of 
V^E'/i follows from Theorem 12.41 provided that we have for all t > 0, 



for any g G C°(M'^) such that /(I + \x\'^)g{x)dx < +oo. Note first that from 
Lemma 13.21 we can restrict ourselves here to test functions g that are com- 
pactly supported. Then we show that the sequence ph is relatively compact 
in C([0, T], LP(R°') — w). This is done by the following lemma: 

Lemma 3.3 For all T > 0, for all p with 1 < p < oo the sequence fh (resp. 
Ph) satisfies 

• fh (resp. Ph) is a bounded sequence in L°°([0, T]; L*'(R'^ x R*^)) (resp. 
m L°°([0,r];LP(R'^)), 

• dtfh (resp. dtPh) is a bounded sequence in L°°([0,T]; W~^''p(U.'^ x R*^))), 
(resp. zn L°°([0,T]; Vr-i'P(R'^)), 

and one can extract from fh (resp. from ph) a subsequence converging in 
C([0,T],LP(R"' X R*^) -u;) (resp. zn C([0, T], Lp(R^) - w)/ 

Proof: the first point uses equation (j^Uj) and Lemma f3. 21 The second point 
uses equation and the identity: 



with the fact that the fh are uniformly compactly supported in x and ^ 
fLemma 13. 2^ : the last point is a classical result of functional analysis (see 



This lemma and Lemma f3. 21 vield ()5Hl . Then using Theorem 12. 4| with p the 
limit of a subsequence of ph, we have convergence of the sequence V'^h to 
V\E'[fi, p] in C([0, T], Lp(R'^)). We have extracted a subsequence fh such that 

• fh converges in C([0, T],Lp(R^ x R"') — w) for every 1 < p < oo. 



(51) 





[T7j for example). 



□ 
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• ph converges in C([0, T], UiW^) — w) for every 1 < p < oo. 

• V\l//i(t, ■) converges in WiW^) for every t and for every 1 <p < oo. 

Thus the hmit (/, V\l/) satisfies equations ()43II44|) and the first part of The- 
orem |S3 is proved. 

3.2 Conservation of energy 

We now give a rigorous proof of the conservation of energy following an argu- 
ment going back to F. Otto (in an unpublished work on the semi-geostrophic 
equations). We recall the definition of the energy as 

E(t) = \ [ fit,x,0\^\^dxd^ + ^ [ p{t,x)\V<i/{t,x)-x\^dx. 

We call the first term the kinetic energy Ec and the second term, multiplied 
by e^, the (normalized) potential energy E^,. We have 

Proposition 3.4 Let f be any solution of \43[ ) such that on every interval 
[0,T], f{t, ■, ■) is uniformly compactly supported in \x\, \^\ < R{T) for some 
function R{T). Then the energy of the solution f is conserved. 

Proof: From Theorem 12.51 we know that 

Ep{t)= I {\x\y2-^>{t,x))p{t,x)dx+ [ {\y\y2-^t,y))dy 



sup / {\x\y2-u{x)) p{t,x)dx+ / {\y\y2 - v{y)) dy 



n 



where the supremum if performed over all pairs {u, v) of continuous functions 
on M'^ such that u{x) + v{y) > x ■ y pointwise. Thus for each t, to G M-I-, we 
have 

Epit) > j {\x\''/2~^>{to,x))p{t,x)dx + j^{\y\^/2-^to,y))dy, 
and this implies 

Ep{t) - Ep(to) > / (|x| V2 - ^(to, x)) (pit, x) - p(to, x))dx 

dtp{s, x) {\x\'^/2 — '^{to, x)) dxds 
e/(s, X, (x - V^(to, x)) dxdids. 



to Jm.'^ 
t 



to JR'ixl 
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Notice that the product in the second hne is hcit since dtp is in W~^'^ for any 
1 < p < oo, fit, ■, ■) and therefore pit, ■) are compactly supported in space 
uniformly on [0, T], and \E' — |xp/2 is in Wl^"^ . Exchanging to and t we would 



have found 

"to 



Ep{to) - Ep{t) > f f ^f{s,x,^){x -V'^{t,x))dxd^ds, 
moreover we have for the kinetic energy 

e^{E,{t)-EM)= [ [ ^f{t,x,0-{^'^{s,x)-x)dxd^ds. 



'to 

Dividing hj t — to,t > to we find 

^, E{t)-Eito) 
t-to 
1 



and 



^ / / ^/(^' ^' ■ (V^(s, x) - V^(to, x))dxd(ds 

t "~ ^0 J to 



^,E{t)-E{to) 



t-to 

< / / ^fis, X, ■ (V^(t, x) - V^(s, x))dxdCds. 

t — to jRd^^d 

We know from l3.ll that V\I'(t, .) converges strongly in Lf^^(]R'^), 1 < p < oo to 
V\E'(to, •) as t goes to to, and so the right hand sides of the above inequalities 
converges to and we conclude that 

t>to t — to 

We could take t < to and find the same result. Finally we conclude that 

dE 



, -0. 
dt 



□ 



3.3 Renormalized solutions and existence of character- 
istics 

The study of renormalized solutions for transport equations has been intro- 
duced in [14] for vector fields in W^'^ with bounded divergence. These results 
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have been extended by Bouchut ^ to the case of Vlasov-type equations with 
acceleration field in BV (A recent result of L. Ambrosio, P], has extended 
the existence of renormalized solutions to transport equations with vector 
fields in BV and with bounded divergence). The fact that solutions of (j43t 
E4|l are renormalized solutions is an immediate consequence of the following 
theorem: 

Theorem 3.5 (F. Bouchut) 

Let f e L°°(]0,T[,L^^(R"' x M"')) satisfy 

§^ + V.-(e/) + V5-(i?(t,x)/)=0, 

with E{t,x) G L\]0,T[;LUR'')) n L\]0,T[; BViU^'')) , 
then, for any g G C^(M), 

^ + V. ■ i^gif)) + • m, ^)9{f)) = 0, 
and for every l<p<oo, f belongs to C(]0, T[, Lf^^(M'^ x M'^)). 

In our case the BV bound on the acceleration V\l' is a direct consequence 
of the fact that ^ is a globally Lipschitz convex function. This result implies 
the strong time continuity results for / and p in Theorem 13.11 Finally, as in 
|14j . it can be deduced from the renormalization property that 

1) for almost every initial condition [x,^), there is a unique trajectory solving 

2) t — > f(t) is just /° pushed forward along these trajectories. 
A complete proof is given in appendix. 

Remark: From the renormalization property it follows that, once the 
potential x) is known, there exists a unique solution to (PHj) in L'^^,^- 
Of course, this does not imply at all the uniqueness of weak solutions to the 
Vlasov-Monge- Ampere system! This paragraph ends the proof of Theorem 

o 



4 Strong solutions 

In this section we show existence of strong solutions over a finite time inter- 
vall. To do so, we need regularity estimates for solutions of Monge-Ampere 
equation. We will get rid of the difficulties that may arise at the free bound- 
ary of the set {p > 0} by considering the periodic case. Note that for the 
Vlasov-Poisson system, existence of global smooth solutions has been proved 
(see (23)) the present case, due to the non-linearity of the Monge-Ampere 
equation, we were only able to obtain a result for finite time. 
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4.1 The periodic Vlasov-Monge- Ampere system 

Polar factorization of maps in a periodic domain 

The polar decomposition Theorem has been generalized in [21] to general 
Riemannian manifolds, while the particular case of the fiat torus T*^ = M'^/Z'^ 
had been addressed in [TT]. 

Definition 4.1 We say that a mapping Y : M.'^ is additive if the 

mapping x —>■ Y{x) — x islf" periodic. The set of all measurable additive 
mappings is denoted by V . For each x we callx the class of x in R°'/Z'^, 
and for any X X the mapping of T"^ into itself defined by 

yx G M^X(x) = Xlx). 

We may say if no confusion is possible additive instead of Z*^ additive. Then 
the following theorem can be deduced from the results of JT] and [21] : 

Theorem 4.2 Let X : M'^ — > M'^ 6e additive and assume that px = X^dx 
has a density in -^^^([0, l]''). Then there exists a unique pair (V$x,7rx) such 
that 

X = V$x O TTx 

where $x is a convex function and ^x{x) — \x\'^/2 isU^ periodic, i^x '-^^ ^ 
^ is additive and %x is Lebesgue measure preserving in T'^. Moreover we 
have 

\\^ — 7rx||L2([o,i]d) = \\X — 7rx||L2(Td) 
and, "^x denoting the Legendre transform of^x, we have 

TTx = V^x o X. 

Remark : The pair ($X7 ^ x) is uniquely defined by the density px = Xjj^dx. 

Notice that the periodicity of ^x{x) — |a^P/2 implies that V$x and V^E'x 
are Z'^ additive, and that "^x — |3;p/2 is also Z'^ periodic. As in the previous 
case, we introduce the following notation: 

Definition 4.3 Let p be a probability measure on T'^, with density in L^{T'^). 
We denote (resp. "^[p]) the unique up to a constant convex function such 
that 

(52) <l>[p] - |x|V2 ^s Z"' periodic , 

(53) yfeC\r), [ fivMp]ix))dx= [ fdp 

(resp. its Legendre fenchel transform). 
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\[^[p] will thus be a generalized solution of the Monge- Ampere equation 
det = p. 

Next the results of Caffarelli (jH], jH], [IDj) on the regularity of solutions of 
the Monge- Ampere equation yield the following theorem: 

Theorem 4.4 Let p > be a C°(T^) probability density on T"^, for some 
a g]0, 1[. 

Then \1' = \l'[p] (see Definition s.^ is a classical solution of 

det = p 

and satisfies: 

II V^(x) - xIIloc < C{d) = Vd/2 
ll^^^llc" < K{m,M, IIpIIc") 

where m = inf p anc? M = sup p. 

This theorem is an adaptation of the regularity results stated above, whose 
complete proof is given in appendix. 



The periodic Vlasov-Monge- Ampere system 

We now seek / : e (T^ x M'^ x [0,T]) f{t,x,0 G M+, for some 

T > 0, solution of the initial value problem for the periodic Vlasov-Monge- 
Ampere {VMAp) system 

(54) % + ^^- (^^) + ■ ((V^M(^) - ^)/) = 

(55) /(0,-,-) = r, 

for a given f^ satisfying the compatibility condition 

(56) j f%x,Odxd^ = l. 

The macroscopic density p is still related to / by equation Q, and \E'[p] is 
as in Definition 14. 31 



4.2 Existence of local strong solutions 

We mention first that the proof of existence of global weak solutions adapts 
with minor changes to the periodic case, and that the obtain for the periodic 
[VMAp) system the same result as Theorem 13.11 
Our result in this section is the following: 
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Theorem 4.5 Let /o G iyi'°°(T"' x M"'), &e swc/i that: 

(57) 3Co > : /o = /or \^\ > Co, 

(58) 3m > : po{x) = / /o(x, 0^^^ > m Va; G T'^, 



then there exists T > and a solution f to the VMAp system (|5^p5,5j) . m 
the space W^''^{[0,T] x T'^ x M^). 

Proof of Theorem 14. 5t First we deduce from Theorem 14.41 
Corollary 4.6 Let p, \1/ = ^[p] be as in Theorem \4.4\ Then, we have 
P'^IIl-CT'^) < C(m,M, llVplUoo(Trd)), 

and we can define 

K{m,M,l) = snp{\\D^^p]\\L^^j,y, ||Vp|Uoo(Tr.) < /, m < p < M} < oo. 

We see that in order to use Theorem 14.41 we need p to be bounded away 
from 0. In the following lemma, we show that under suitable assumptions 
on the initial data, it is possible to enforce locally in time the condition 
<m < p < M. 

Lemma 4.7 Let f G L~([0,T] x T'^ x R"^) satisfy 

^ + V.-(e/) + V5-(i5;(t,x)/) = 

/(o,.,.) = r, 

with E G L\[0,T],BV{T'^)) and 

ll-^llL°°([o,r]xT'*) < 
let the initial condition /o be such that 

a(x,O</(0,a;,O<Ka:,O, 

with Pa{x) = J a{x,^)dC, > m > and Pb{x) = Jb{x,^)d^ < M < oo and 
a, b satisfying 



c 



|V.,«(a,6)|<^^ 

Then there exists a constant R > depending on m, M, c, F, such that 

(Paix) - Rt) < p(t,x) < (pb(x) + Rt). 
The proof of the lemma is given in appendix. 
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4.2.1 Construction of approximate solutions 

Let us consider {t, x) — * E{t, x) a smooth vector-field on T*^, and write 

If / satisfies TE{f) = 0, we have 

TEdtf = -dtE ■ VJ, 

and therefore 

(59) ^||V.,5/|Uoo < ||V.,5/||loo(1 + 
which implies 

||V,,5/(t)|Uoc < ||V,,^/(t = 0)|Uoo exp (^^ (1 + \\V,E{s)\\Lo^)ds^ . 

Now let /o be given as in Theorem 14.51 satisfying ()57l58p . Thanks to Lemma 
14.71 it is possible to find ti,m,M such that for any / satisfying 

TEif) = 
fit = 0) = /o, 

for any field E G L\[0,ti], BV{T'^)) satisfying ||E||ioo([o,t,]xT<') < Vd/{2e^), 
we have 

(60) m < p(t, ■) < M, Vt G [0,ti] 

(61) leUax < Ci = lOCo, 

with / supported in {\^\ < \^\max} and with Cq as in Theorem 14.51 so that 
we have for < t < ti: 

(62) llVplUoo <a;dCf||V./|Uoo, 

Ud being the volume of the unit ball of M''. Then we construct a family of 
approximate solutions {fh, "^h) to (|^. in the same spirit as we did for weak 
solutions, by solving 

f + e-v.A + ^^^^-v,/.^o 

fh{t = 0) = /o 

\[/^(t) = \[/[p(t = nh)] for t G [nh, {n + l)/i[. 
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Note that we have neither molhfied the term V^E' nor the initial condition 
and that || V\l//i — xH^oo < C{d) = \/d/2. Now choose I = lO\\Vx,(fo\\L°°'^dCf, 
if for some t = nh < ti — h we have 



|V.,g/"(t = n/i)|Uo. < ^ 



this imphes, thanks to (j62|) . that 

\\V,p\t = nh)\\Lo. <l, 

and conditions ()60p61|) are satisfied because t < ti. Then if we denote K = 
K{m, M, I) as in CoroUarv 14 ■6[ we have for nh < t < nh + h, 

and then 

= nh + < W^.^i fit = nh + h)\\L^ exp (K + l)h. 

So if we define T as 

T = min{ti,t2}, 
with exp((i^' + l)t2) = 10, we have for < t < T, 

II Vx.,^/^||loo < 10|| Va:,^/o||L°° 

IIVp'^iU- < / 

m< p< M 
P^^^IU- < K. 

Thus we can extract a subsequence converging to a strong solution of 



()54p55|l . Then we argue as in section 2 to show that all terms converge to 
the correct limit. This ends the proof of Theorem 14. 51 

□ 



5 Asymptotic analysis 

5.1 Convergence to the Euler equation 

In this section we justify that the Vlasov-Monge- Ampere system describes 
approximate geodesies on volume preserving transformations: indeed we will 
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show that weak solutions of this system converge to a solution of the in- 
compressible Euler equations (E) as the parameter e goes to 0, at least for 
well prepared initial data. We restrict ourselves to the space periodic case, 
the macroscopic density p is still defined by (jH)) and the convex potentials 
\E'[p] are still as in Definition 14.31 
For sake of simplicity, we slightly modify our notations and introduce the 
following rescaled potentials 

|x|V2-^[p] 



e 

|2 



KM - |x|V2 

np\ = ^ , 

so that 

V<^[p] = V<^[p]oV$[p], 
and the VMAp system takes the following form: 

(63) | + ^.v./-Md.v,/ = 

(64) /(O, ■, ■) = /„. 
The energy is given by 

(65) E{t) =^ I f(t,x,0\^\'dxd^+^ I iW^l'dx 



I J fit,x,m'dxd^+l J plV^prfx. 



It has been shown in section 13.21 that the energy is conserved. The Euler 
equations for incompressible fluids (E) reads: 

(66) dfV + V ■ Vv = —Vp 

(67) V ■ = 0. 

We shall here consider a smooth solution of E and a weak solution of VMAp, 
with 'well prepared initial data', meaning that the initial data of both systems 
are close a time 0. Then we will show that as time evolves, both solutions 
stay close to each other. 

Theorem 5.1 Let f be a weak solution of iS^A with finite energy, let 
(t, x) v{t, x) be a smooth C^{[0, T] x T'^) solution of ^E3^ for t G [0, T], 
and p{t, x) the corresponding pressure, let 

He{t) = \ I f{t,X,i)\i~v{t,x)\^dxdi+]^ I \Vv{t,x)\^dx, 
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then 

H,{t) < Cexp(Ct)(/f,(0) + e^), Vt G [0,T]. 
C depends only on T,s\x^Q^^^r^ {\\v{s, .),p{s, .),dtp{s, .),Vp{s, .)\\who.(jd)} . 

Remark 1: This estimate is enough to compare the weak solutions / to 
the VMAp system (for well prepared initial data) and the smooth solu- 
tions V of the Euler equations. For instance, / f{t = 0,x,^)d^ = 1 implies 
ip{t = 0,x) = and therefore, 

y le - v{t = 0, x)\^f{t = 0, X, Odxd^ < Coe' 

implies 

sup / \^-v{t,x)\'^f{t,x,^)dxd^<CTe'^, 
te[o,T] J 

where Ct depends only on Cq, T and v. 

Remark 2: We see that we consider nearly monokinetic initial data for the 
VMAp system. 

Proof of Theorem 15.11 



H.= - j fit,x,m-v)\/vi^-v) 
+ j f{t,x,i)-v-V^ 



We shall show that 
d_ 
~dt 



(68) - j /(t,x,0(t^-0-Vp, 

where we will use the notation 

d 

u Vf w = u^diV^w^ . 

The proof of this identity is postponed to the end of the section. 

Now we look at all terms of the right hand side. All the constants that 
we denote by C are controlled as in Theorem 15.11 We set 



T, = - [ /(t,a;,0(^-0-Vp. 
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First we have Ti < CH^. For T2 we have 



To 



pv ■ Vip 



t;(V$[p]) ■ V(^(V$[p]) 



1 f 

= - V ■ Vip + {v{x + eVifi) — v{x)) ■ Vip 



we have used that v is divergence-free thus its integral against any gradient 
is zero. Next we have the following lemma: 

Lemma 5.2 Let G : T"' ^ M &e Lipschitz continuous such that / G = 0, 
then for all R > 0, one has 

pG\ < ^\\VG\\Loo{^e' + RH,). 

Proof: We just write a Taylor expansion of G: 



(p-l)G 



{G{x + eV(^) - G{x] 



< e\\VG\\L^\\Vip\\Li < e\\VG\\L^Hl/^ < h\VG\\L^{^e^ + RH, 



□ 



Again, since v is divergence-free, J v ■ Vp = 0, thus from Lemma [5.21 we have 

- ^ pvVp<Gie^ + H,). 



We remind that 



dtpit,x) = -V,.- j 



Since it costs no generality to suppose that for all t G [0,T], J p{t,x)dx = 0, 
we obtain that 



/(t,x,oe-Vp 



dp 

d f [dp 

= Jt]""'j"Tt 
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again using Lemma f5. 21 where Q{t) = — J pp. Thus 



2^ dQ 



dt 

and we have the following inequality: 

(69) ^^(H, + Q)<CH, + 0{e'). 

Moreover, using Lemma f5 .21 

(70) |Q(t)| <Ce2 + /7,(t)/2, 
thus 

(71) H, + Q> HJ2 - Ce\ 
and we can transform ()69|) in 

(72) j^{H, + Q)< C{H, + Q) + Ce\ 

Gronwall's lemma then yields 

HA^) + Q{t) < iH,{0) + g(0) + Cte^) exp(Ct). 
Using again (f70|) we obtain 

(73) H,{t) < CiH,{0) + e') exp(Ct), 
which achieves the proof of Theorem 15.11 



Proof of identity 

We first notice that, for all G C-^(R x T"^) , we have: 

^ J pit,x)g{t,x)dx = J J f{t,x,C){dtg{t,x) + i ■yg{t,x))didx. 
We also use the conservation of energy defined by (|65|l . Then we get 

jH.= ~ I fit,x,mv\'-2^-v)dxd^ 

= j f{t,x,0{dtv-v-dtv-i)-]^ j V. ■(/(t,x, 00(1^1' 

+1 j v^.{Uj^f{t,x,mv?-'2i-v). 

25 



Integrating by part, we get 

d 
di 



The first two terms can be rewritten as 

J f(t,x,0{dtv-v-dtv-0 + J /(t,x,OeVt;(t;-0 

f(t,x,Oiv-0^viv-0 + j fit,x,OdtV-{v-0 
+ J fit,x,Ov^viv~0 
= -j fit,x,Oiv-mviv-0 + J fit,x,Oiv-0-idtV + vVv), 
and finally using equation (j66|l we conclude. 

□ 

5.2 Comparison with the Euler-Poisson system 

Here we show that, as mentioned in the introduction, the Euler-Poisson (EP) 
system is a more accurate approximation to the Vlasov Monge-Ampere sys- 
tem than the Euler equations, as e goes to zero. 

The EP system Let us recall that the (pressureless) Euler-Poisson system 
describes the motion of a continuum of electrons on a neutralizing background 
of ions through electrostatic interaction. Let v and p be the velocity and 
density of electrons. Let (f be the (rescaled) electric potential. Under proper 
scaling, these functions of x G M'^ and t > satisfy the Euler-Poisson system: 

(74) dtv + v - Vv = --V(p 

e 

(75) dtp + V- (pv) = 

(76) 1 - eA^ = p. 

The so-called 'quasi-neutral' limit e ^ of similar systems has been studied 
for example in ^H] and jT^j, and convergence results have been established 
using pseudo-differentials energy estimates. For well-prepared initial data, 
solutions of EP are expected to behave as solutions of Euler incompressible 
equations. This fact is proved by the second author in his PhD thesis (|l2n]. 
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Chap 2), see also [1^]. We give here the complete result that we will use 
herafter. We will denote by v"^ (resp. f^) the solutions of the EP (resp. 
VMAp) system with parameter e. 

Theorem 5.3 Let v be a solution of ^6(^61}) on [0, T] x T^, with initial data 
Vq, and satisfying v G L°°{[0, T], if^(T'^)) for some s > So{d). There for some 
s' > 0, s' < s, if {vq,Pq) is such that the sequences 

Vq-vq pi - 1 

e e-^ 

are hounded in H^'iT'^), then there exists > with liminfe_^o^e > T and 
a sequence {v'^,p'^) of solutions to the EP system on [0,T^[ with initial data 
{vq,Pq), belonging to L°°{[0,T^], H'"^ i'^'^))- Moreover, fore small enough, the 
sequences 

v"^ — V ~ 1 

are hounded in L°°([0, T], //^'(T'^)). Finally, s' goes to +oo as s goes to +oo. 



Assumptions 

Here we consider v a solution to E (IMl IHZj) with initial data fo, a sequence 
f of solutions of VMAp ()63|64|1 with initial data /q, and a sequence {v^.p") 
solutions of EP (fflj IT5| I76|l with initial data (vg^Po)- We still define as 
in Theorem 15.11 

Heit) = \ j r{t,X,m-v{t,x)\''dxdi+]^ j \Vip'\^dx. 

We introduce the following assumptions: 

HO V solution of E satisfies, for some Co > 0, llfH^ooQQ j^jj^^scj-d)) < Co, and 
s is large enough so that s' in Theorem 15.31 satisfies s > s' > [^] + 2. 

HI The sequence {vq, Pq) of initial data of EP is such that, for some Ci > 0, 

sup I -\\Vo-v\\jjs'„d), \\\p' - l||//.'(Td) \ < Ci. 

H2 The sequence /q satisfies H^{0) < €2^^ for some C2 > 0. 



HO, HI, H2 imply that 



27 



1. There exists Cq such that 

(77) II <Co. 

2. From Theorem 15. 11 there exists C*i such that 

(78) H,{t) < Cie^ ioT t e [0,T]. 

3. From Theorem 15.41 and Sobolev imbeddings, there exists C2 such that 
v'^ — V P*^ — 1 



(79) 



sup 

e<eo 



< Co 



L°°([0,T],H/2,oo(T<*)) 



We are now ready to prove the following result: 

Theorem 5.4 Let f^yV^, Pq,v,T be as above, satisfying assumptions HO, 
HI, H2. Define 

Ge{t) = l J nt,x,0\^-v'{x)\'dxd^ + \ J iVy^^-V^fcix. 

Then there exists C > such that 

G,{t) < Cexp{Ct){G,{0) + e^), Wt e [0,T] 
where C depends on s', Co, Ci, C2, T. 



Remark: the theorem shows that the distance between solutions of the {EP) 
system and the VMAp system measured with is like O(e^) whereas Theo- 
rem showed that the distance between the solution of the Euler equation 
and the VMAp system was like O(e^). Note also that G^ and can both 
be interpreted as the square of a distance. 

Proof of Theorem 15. 4t For notational simplicity, we drop most e's. Pro- 
ceeding as in (j68|l and noticing that: 

4 / |V<^|' = - / pv-V^ 
ax jjd e jjd 
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we obtain the following identity: 
d 



^0) -T" / V<^ ■ V<^. 



dt 



Then we notice 



fit, a;, 0-^ ■ V<^ = ^ ( /" -pip] - - I pdtif. 



e dt \J e J e 

Next, we have the following lemma: 
Lemma 5.5 Define for any 9 E C^(T'^) 



<V9 > (x) = I Ve{x + seVip{x))ds, 
Jo 

<V'^0>{x)= [ {l-s)V'^e{x + seVip{x))ds. 







Then 

^ p9 = I 9 + e I <\/9 > -Vifi 



9 + e j V9-V^ + e' j <V'9> V<^V<^. 

Proof: The proof just uses the Taylor expansion and the identity 
^ p9 = ^9{x + eW^). 

Using Lemma 1331 we get 
j pdt(p 

We claim that, under our assumptions, we have 

ll"9tV^(^||Loo([o,T']xTd) < C. 
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□ 



Proof: from (fTSj) . we have 

dtp = —pv ■ V — V ■ Vp. 

Using (f?^. we obtain that ||5tp||j:^s'~i < Ct. Since H'^' {T'^) is continuously 
embedded in W'^'°^{T'^), H''-\T'^) is continuously embedded in L°°(T'^). 
Then, using (fTBj) and classical elliptic regularity, we have 

eWdtW'^ipWHs'-i < C\\dtp\\Hs'-i, 
and the desired result follows. 



□ 



This implies, using fITSj) . that 



e / < dtV'^ip > VifVf 



Next, 



Using again Lemma 1^31 we get 



e dt 

and for the same reasons we have || V^<^||Loo([o,T]xTd)) ^ Ce. This yields 

Q{t) = t j <V^^> VipVifi = O(e^). 
Moreover, it does not cost to set J ip = and deduce 

I f{t,x,0k-^V-j^ |v^-V^ = -^y" pv.Vv + 0{e') + j^Q. 
Thus the remaining terms are 

R = - J [pV(^ — pV(^ + pV(^ — pVv?] ■ V. 
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Calculations that we postpone to the end of the proof show that 
(81) ~^ / ^ ' ^^'^^'^ ~ ■ '^^^ + 

with C depending on || V^-uHlooqo^tIxT'^) and e~^|| V^<^||2,cx)([o,T]xTd), therefore 
uniformly bounded thanks to (f79|) . Finally we obtain 

jG,< - j /(t,x,0(e-^)Vt;(e-t)) + (V^-Vv9)V^;(V^- V<^) 
/(V ■ v){\^^? - 2V^ ■V^) + C j - V^p 

with |g(t)| < Ce'^ for t G [0,r]. From ^ we have ||V ■ t;||Loo([o,T]xTd) < C'e 
and II V<^||j;^oo([o,T]xT'*) ^ C'e, whereas (|78|l yields y |Vv9|^ < Ce^. Note that 
we also have 



- j fit, X, 0(e - v)^v{^ -v) + (V^ - V^)Vv{V^ - V<^) 
+C / |V<^ - V<^|' < CG,. 



We conclude that 

|(G.-g)<C((G,-g) + e3), 
and the conclusion of Theorem 15.41 follows by Gronwall's lemma. 

□ 

Proof of identity (j81|) : Here we have to compute: 

R = - I v{x + eVif) ■ Vip - {vV^){x + eVif) + (1 - e/\Lp){v -Vip-v- Vip) 



e 

Using Lemma f5. 51 we have: 
If 

R = - V ■ Vv9 — V ■ Vy? + V ■ \/(p — V ■ Vy? 
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We see that the first hne cancels. Then we show that the last line is bounded 
by Ce\ 

This is obvious for the last term since from (|77| I78|) we have H-uH vi/2>°° < C, 
and II V<^|| vi/2,=o < Ce. 

Then for the first term we have the following lemma: 
Lemma 5.6 We define 

then one has: 

|A| < Ce^°/3 + C J |V^ - V^l^dx. 
Proof: First we show that if Q{R) = /{|v^|>_r} I^V'I^; 

e{R) <c j |v^- v^|^ + 

Proof: J iVv^l^ < C'e^, implies that 



meas{|V^| > R} < C{^f. 

R 



Since \V^{t,x)\ < e for {t,x) G [0,T'a;] x T'^ 

Q{R)< [ |V^|'+ / \Wip-Vip\ 

J^\V^\>K\ J\\Vio\>R\ 



4 



2 



Then we have 



A < Ce{R) + [ \<Vv> (x) - Vv{x)\ VipVip 



|V<p|<i? 

with |< Vt; > (a;) - Vt;(a;)| < Ce|V(^| 
thus A<Ce [ + Ce{R) 



2 , ^ , / n-v ^ -\2 



<C[eR I \Vif\' + — + I iVv^-V^I 



<C[e'R + +—+ I iVy.- V^r 
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for all R, so for R = e^^^^^ one obtains: 



□ 



Thus we have shown that R — S + O(e^), and S — E^^^Tfe where each is 
given by: 

Ti = djVidj(pdi(p 
T2 = -djVidj(pdi(f 
T3 = -Vidijipdjif 
T4 — djVidj(fdi(f 
T5 = Vidij(pdj(f 
Te = Vidjj(fdi(p 



where we have used Einstein's convention for repeated indices. First we have 

Ti+T2 + T4^ J djViidjif - dj<f){di(f - di(p) + Ty 
with Ty — J djVidjipdi(f. 



and 



thus 



Te = J -{V ■v)A^^ + Ts-T3 
with Ts — J djVidijifxp. Then 

Ts= - djVidj(pdi(p + dijVidjip(p 



-T7 + / V ■ v{Aifif + WifWif) 
33 



and finally we obtain 

S{t) = I VD(V<^ - V<^)(V<^ - V^) - ^(V • v) |V<^ - Vy?!^ 



and the identity (j8T| is proved. 



□ 
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6 Appendix 



Existence and uniqueness of solutions to second order 
ODE'S with BV field 

In this section we prove existence and a.e. uniqueness for ordinary differential 
equations of the form: 

for X e T'^, Ye T'^, and where the field E belongs to L°°(]0, r[xT'^) n 
L^(]0,T[, BV(T'^)). We work in the flat torus for simplicity, but our results 
are still valid in an open subset of W^. This result is an adaptation of the 
proof of that uses the result of |0] on renormalized solutions of transport 
equations. 

Remark: After this proof was written, the authors learned of a result by L. 
Ambrosio ([Tj) that extends the results of [Hj to transport equations when 
the vector field is in BV with bounded divergence. 

Renormalized solutions for Vlasov equations with BV field 

Theorem 3.4 in |^ adapted to the periodic case sates that if / G L°°{]0, T[xT'^x 
R"') satisfies: 

(83) ^ + ■ + ^« ■ = °' 

with E{t,x) e Li(]0,T[xT"') n L^(]0,T[, BV{T'^)), then for all g Lipschitz 
continuous we have 

^ + V. ■ mf)) + ■ m, ^)9{f)) = 0. 

The property of renormalization implies that 

• solutions to (j83|l with initial data in ^^^(T'^ X M'^) belong to 
C(]0, T[, Lf^,(T^ X R"')) for any 1 < p < oo, 

• solutions to ()83p with prescribed initial data in L°°[T'^ x M*^) are a.e. 
unique. 
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• if En converges to E in L^(]0,T[xT^) then the solutions of (jHH)) with 
En instead of E converge to the solution of (jH^ . 

We notice that equation (j54|l satisfies the assumptions of the Theorem, and 
thus will have the renormalization property. This renormalization property 
was used in jll] to obtain a.e. uniqueness for solutions of the corresponding 
ODE'S. Indeed, the ODE 

dtX{t,s,x) = b{t,X) 
X{s, s,x) = X 

is associated to the transport equation: 

dtu + b{t,x).Vu = 

whose solutions satisfy for all (t, s) g]0,T[ 

u(t, X(t, s, x)) = u{s, x). 

We extend this consequence to the case of second order equations, with BV 
acceleration field. To the kinetic equation 

(84) oj + ^.^j + E{t,x)-VJ = 

we associate the second order ODE (j82|) which can rewritten as duX = 
E{t,X). The result is then the following: 

Theorem 6.1 Let E{t, x) e L~(]0, T[xT^) n L\]0, T[, BV{T'^)), 
then the ODE 

(85) duX{t,s,x,0=E{t,X) 

(86) (X(s, s, X, 0, dtXis, s, x, 0) = (x, 

admits an a.e. unique solution. 

Remark: Here almost everywhere must be understood for the Lebesgue mea- 
sure of M^. 

Proof of Theorem l6.lt We know that through equation (j82|) equation (j85|l 
can be considered as a first order differential equation. Let us first consider 
the case where E is smooth. Note y G T"^ x M'^ (resp. y) for {X, V) (resp. 
for (a;,0) and B eR'^xR'^ for i^,E). Then for all s g]0,T[, Y solves: 

(87) dtYit,s,y) = Bit,Y{t,s,y)) 

(88) Y{s,s,y) = y 
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Then for all t, ti, ^3 g]0, T[ we have the following: 



Y{t,t,y) = y 
Y{h,h,Y{t2,h,y)) = y. 



Differentiating the last equation with respect to ^2 yields: 



(89) 
(90) 



dsY{t,s,y) + VyY{t,s,y)-B{s,y) = 
Y{t,t,y)=y. 



Yt{s,y) = Y(t,s,y) thus solves a transport equation which is nothing but 



Now we show existence and uniqueness for solutions of ()87|88|) . Let t and 
s be fixed. Let us consider a regularization En of the the field E and set 
Bn= (e,^n). We have 

• t Yi,„(t, s,y) that satisfies (|87I88|) 

• F2,n(t, s,y) that satisfies (j89l9n|) . 

From the stability Theorem 2.4 in we know that the whole sequence 
t F2,n(t,s,.) converges in C(]0, T[, Lf„^(M'^ x T"^)) to t l"2(t,s, .), the 
unique renormalized solution of ()89I90|) . Thus for fixed t the whole sequence 
Y2^n{t,s,.) converges strongly in Lf^^{M.'^ x T"^). Now since for every n we 
have Yi^n(t,s,y) = Y2^n(t,s,y) the same property holds for Yi^„(s,t, .). Now 
we can pass to the limit in the term Bn(t, Yi^nit, s, y)). Indeed, by density of 
functions in L^, if we have Es G approximating E then 



The second term goes to because of the strong convergence of Yn, the first 
and the third go to because Y and Yn are measure preserving mappings, and 




\\B{t,Yn{t,s,y))-B{t,Y{t,s,ymLi 
< \\B{t, r„(t, s, y)) - B,{t, Yn{t, s, y)) lUi 
+ \\Bsit, F„(t, s, y)) - Bs{t, F(t, s, y)) |Ui 
+ \\B{t,Yit,s,y)) - B,{t,Yit,s,ymLi 
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so for example \\B{t, Y{t, s, y)) - B,{t, Y{t, s, y)) |Ui = \\B{t, y) - B,{t, y) l^i. 
So finally we have 

\\B,,{t,Yn{t,s,y)) - B{t,Y{t,s,y))\\Li 
< ||5„(t, s, y)) - Bit, F„(t, s, y))\\Li 
+ \\B{t,Y^{t,s,y)) - Bit,Y{t,s,ymL^ 

that goes to and we can pass to the limit in the equation (j87p88j) and the 
existence of a solution to ()87p88|l is proved. 

To obtain uniqueness, we argue as in [14 . Any function of the form go(Y(t, s, y)) 
is a solution of ()91I92|) . thus by uniqueness of the solution of the transport 
equation we obtain uniqueness of the ODE. 

□ 

A remark on ODE's of second order 

In this section, we want to solve the Cauchy problem for: 

dttX{t,x)=E{t,X) 
(X(0,x),9iX(0,x)) = (a;,t;(a;)) 

with E as above. We are thus interested in monokinetic initial data. 

Theorem 6.2 for all v^{x) vector field on T'^ , and for Lehesgue almost every 
Sv G M.'^, there exists an a.e. unique solution to 

duX{t,x) = E{t,X{t,x)) 

(X(0, x) , dtX{0, x)) = (x, v^{x) + Sv) 

Proof: Let g{x,^) be the indicator function of the set of those [x]^) such 
that the trajectory coming from x is not well defined. We just have to prove 
that for a.e. 6v G M'^ we have J g{x, f "^(x) + Sv)dx = 0, which is true because 

g{x, v°{x) + i)dxdE, = j g{x, i)dxdi = 0. 

Stability Using the fact that for En converging to E in with 
E G L\]0, T[, 5\/(T^)), we have X„(t, x, v) X(t, x, v) in C([0, T], L^), we 
have then, for all t, for almost every 6v, X„(t, x, v^{x) + 6v) X{t, x, + 
6v) in L^. Thus we have 
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Theorem 6.3 If E.^ converges to E in let he solution of 

dttXnit.x) = En{t,Xn{t,x)) 

(X„(0, x), 9tX„(0, x)) = (x, v\x) + 5v) 

then for all t, for almost every 5v, X„ converges in L^(R'^) — s to a solution 
(unique for almost every 6v) of 

duX{t,x) = E{t,X) 

(X(0, x), dtX{0, x)) = (x, v\x) + 6v). 

Control of macroscopic density in kinetic 
equations 

We prove here Lemma [4.71 

Lemma 6.4 Let f G L°°{[0,T] x T'^ x R'^) satisfy 

(93) ^ + V.-m + V^-{E{t,x)f) = 

(94) /(0,.,.) = r 
with E G L\[0,T];BV{i:^)) and 

(95) II-^IIl°°{[o,t]xT'') < 
Let the initial condition /q be such that: 

a(x,O</(0,x,O<Ka:,O, 

with Pa{x) = J a{x,^)d^ > m > and Pb{x) = Jb{x,^)d^ < M < oo and 
a, h satisfying 

Then there exists a constant R> such that 

(Paix) - Rt) < p{t, x) < (p,(x) + Rt). 

Proof: First suppose that the force field and the initial data are smooth. For 
equation ()93I94|) we can exhibit characteristics {x,^){t;tQ,XQ,^o), giving the 
evolution of the particles in the phase space. We have f{t, x, ^) = f{tQ, xq, Co)- 
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Since the initial data is compactly supported and the force field is bounded 
in the norm, we have 

\^-^o\<F\t-to\, 

\x-xo\ < i\U + j\t - to\)\t - to\. 
If for t = we have a(x, ^) < /(0,x, ^) < b{x,^) then 

A{t,x,0 < f{t,x,0 <B{t,x,0 

F 

A{t,x,i)= inf a{x + \t-tQ\{i + -\t-tQ\)ai,i + F\t-U\G2) 
_ F 

B{t,x,i)= sup h{x+\t-tQ\{i + -\t-tQ\)ai,i + F\t-tQ\a2). 

ki|>k2|<i ^ 

Using and integrating in ^ we find thus a constant R = R{F, C, d) such 
that for t — to ^ 1 we have: 

Pa{x) - R\t - tol < p{t,x) < pb{x) + R\t - to\. 

Next we need to show that the solution of the regularized equation converges 
to the solution we are studying: this result comes from the uniqueness of 
the solution to ()93p94|l which is a consequence of the renormalization prop- 
erty. Indeed since E is bounded in BV the system ()93I94|) admits a unique 
renormalized solution and the sequence of approximate solutions converge in 
C([0, T], for 1 < p < cxD thus the bounds obtained above are preserved. 

□ 

Regularity of the polar factorization on the flat torus 

Here we deduce from HH and |H1, 0, [IDl the Theorem lOl 

Theorem 6.5 If p e C"^(T'^) with <m < p < M is a probability measure 
on T'^ then \1' = \l'[p] (see Definition \4.3}) is a classical solution of 

(97) detL>2^ = p 
and satisfies: 

(98) II V^(x) - x|Uoc < C{d) = Vd/2 

(99) jjl^^^llc- < K{m, M, ||p||c«) 
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Proof of Theorem 16.51 Consider p a Z'^ periodic probability measure, 
satisfying 

(100) <m < p < M, 

and $[p] as in Definition 14.31 First it is shown in ^T] that 

(101) \V^p]{x) -x\ <C{d). 

It follows that the strict convexity argument of 1^ applies: indeed if $ = $ [p] 
is not strictly convex its graph contains a line and this contradicts poip . 
Moreover since $ — |xp/2 is globally Lipschitz and periodic there exists N{d) 
such that 11$ — |xp/2||ioo < N{d). It follows then that there exists < 
r{d) < R{d) and M{d) such that 

(102) B{r{d)) C {$ - $(0) < M{d)} C B{R{d)) 

It remains to show that our solution is a solution in the Aleksandrov sense 
of the Monge-Ampere equation 

m < det D^^ < M. 

This is a direct consequence with minor changes (to adapt to the periodic 
case) of Lemma 2 of [10 . Then, normalizing $ to $ = $ — $(0) — M{d) it 
follows that $ is a solution of 

p(V$)detL'^$ = 1 

$ = on on 

B{r{d)) cnc B{R{d)) 

Thus the interior regularity results of pj apply uniformly to all $[p] with p 
satisfying ()100|1 and ||p||c"{T<^) bounded and Theorem 14.41 follows. 

□ 
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